Helioseismology with PICARD 
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Abstract. PICARD is a CNES micro-satellite launched in June 2010 [17] . Its main goal is 
to measure the solar shape, total and spectral irradiance during the ascending phase of the 
activity cycle. The SODISM telescope onboard PICARD also allows us to conduct a program 
for helioseismology in intensity at 535.7 nm [5]. One-minute cadence low-resolution full images 
are available for a so-called medium-/ program, and high-resolution images of the limb recorded 
every 2 minutes are used to study mode amplification near the limb in the perspective of g-mode 
search. First analyses and results from these two programs are presented here. 



1. Medium-/ program in intensity 

Full continuum images at 535.7 nm (bandwidth 0.5 nm) are recorded by SODISM [8] every 
minute with a spatial resolution of about 1 square arcsecond on a 2048 2 CCD. For telemetry 
reasons, these images are downgraded to 256 2 pixels (about 8x8 arcsecond) by a simple 
onboard binning before being transmitted to ground. We have analyzed a dataset covering 209 
consecutive days (2011/04/16-2011/11/10 ) with a duty cycle of 74.4%. The gaps are mostly 
due to interruptions for the other PICARD program (astrometric measurements at different 
wavelengths) and to the crossing of the South Atlantic Anomaly which is the location of transient 
signals generated by trapped protons exposing low-altitude, Sun-synchronous orbiting satellites, 
like PICARD, to strong radiations. 

Figure [I] (top) shows the l-v diagram obtained from this 209-day dataset for angular degree 
up to Z=300. As expected from the instrumental resolution, a spatial aliasing for spherical 
harmonics above Z=200 is present. We can also see a temporal aliasing for frequencies above 4.5 
mHz. This is the signature of a spurious cycle in the 2-min sampling slots corresponding to the 
astrometric program sequences. Figure [I] (bottom) illustrates a m-v diagram for Z=90 clearly 
showing the S-shape signature of the differential rotation. Figure [2] shows the mode parameters 
modelled using the asymmetrical line profile formula of Nigam and Kosovichev [10] and fitted 
using a Maximum Likelihood Estimator minimization [HE]. For a given mode (n,/), the 2/+1 
m-components of the multiplets are fitted simultaneously using an expansion on orthogonal 
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Figure 1. Top: l-v diagram 
from 209-day PICARD observa- 
tions. The fitted modes are shown 
by the green dots. Bottom: vn-v 
diagram for Z=90 and radial orders 
n =7 (left) and n=8 (right). The l- 
leaks from Z=90±3 are clearly seen. 



Figure 2. Fitted mode parameters: amplitude (upper 
left), linewidth (upper right), peak asymmetry (lower 
left), and as splitting coefficients (lower right). The red 
line is a fit of SODISM positive asymmetries in intensity, 
while the green line shows the negative asymmetries 
obtained from the SOHO/GOLF velocity signal. 



polynomials [12] with up to 9 so-called a-coefflcients. The / and m mode leakages are included 
in the fitting process. Positive peak asymmetries are obtained as expected from previous analyses 
of intensity signal. The as splitting coefficients plotted as a function of the mode turning point 
proxy i//L, where L = + 1), exhibit a gradient which is a typical signature of the solar 
tachocline. Improvements of the peak-fitting process are currently underway. We need for 
instance to better model all the components of the leakage matrix in order to fit higher angular 
degrees. 




500 



450 



400 i 



350 



300 



Figure 3. Symmetric part of 
the internal solar rotation rate 
with respect to the equator. Iso- 
contours are shown every 10 nHz 
from 300 nHz up to 480 nHz. 
This was obtained by inverting 
the frequency splittings estimated 
from PICARD/SODISM intensity 
images over the period from April 
16 to November 10 2011. The 
solution below O.3i? is poorly 
constrained by the data and is not 
shown. 



Figure [3] shows a Regularized Least-Square (RLS) inversion [6J of the measured rotational 
splittings of 769 modes from l—l to Z=99. The steep gradient of the tachocline at the base of 
the convection zone around 0.7 Rq is clearly seen. The inferred rotation rate is similar to the 
one previously obtained from velocity data [6j [TT] : a radial gradient is obtained close to the 
surface followed by a latitudinal differential rotation persisting through the whole convection 
zone. The radiative interior is roughly rigid at a rate of 435 ± 5 nHz which corresponds to the 
surface rate for latitudes between 30 and 40 degrees. Finally, we note that a local maxima is 
found close to the equator near O.95i? . These preliminary results are very encouraging. We are 
currently working in order to extend the peak-fitting analysis towards higher / and to increase 
the number of fitted a-coefflcients thus improving the resolution of our inversion. The reversal 
of peak asymmetry observed between velocity and intensity power spectra is linked to the effect 
of a background component correlated to the modes [9] . The complete two years of PIC ARD 
observations available to date should then be used to confront the details of internal rotation 
studies obtained using simultaneous velocity and intensity data from SDO/HMI or GONG. By 
fitting intensity, velocity, their phase difference and coherence signals simultaneously, we hope 
to nail down systematics in either observables and better understand the correlated background 
noise [HJH]. This will help to better interpret the fitted frequencies and frequency splittings in 
terms of internal structure and dynamics. 

2. Limb seismology 

Every 2 minutes, the full resolution of PICARD/SODISM images are kept in an annulus of 
22 pixels around the solar limb (increased to 31 pixels since April 2012). In April 2011 and 
April 2012, 3 consecutive days were dedicated to helioseismic measurements at 535.7 nm without 
interruption for the other PICARD objectives. These datasets with a duty cycle above 95% have 
been used in order to study the expected mode amplification at the extreme limb [21 \W\ [T5] . 
Intensity data from SOHO/MDI and SDO/HMI were also used for comparison. 
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Figure 4. Classical spherical harmonic transform of full images {l-v diagram ) (left) and limb 
transform of limb pixels (I'-v diagram) (right). The theoretical n—Q ridge is indicated by a 
dashed line for both cases clearly showing a frequency shift between the two diagrams. 

In order to analyse the limb signal, we calculated the Fourier transform e ll ' p along the limb 
using the polar coordinates, where p—1 at the disk center and p=0 at the limb. Figure [4] shows 
the comparison between the l-v diagram obtained from a classical spherical harmonic transform 
of full images (left) and the l'-v diagram (right) resulting from the limb transform as described 
above. A frequency shift between ridges of the same radial order is clearly seen. We carried out 
simulations in order to understand this shift and the link between the two diagrams [3J. For a 
given it can be shown that the limb response are non zero for value of / higher than V . Since the 
mode frequencies increase with the degree, the leaks will appear at higher frequencies than the 
original anticipated frequency. The leaks will then appear as an asymmetry at higher frequency 
with a decreasing amplitude when I increases. This is indeed what is shown in Figure [5j where 
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Figure 5. Comparison between observed and simulated profiles. On the left panels, the 
black lines are ridges in SOHO/MDI l-v diagrams while the green lines are ridges in the 
PICARD/SODISM (top) and SOHO/MDI (bottom) V-v diagrams. The two panels on the 
right show the simulated profiles for both the l-v (black) and V-v (green) diagrams which are 
also shown with an appropriate scalling on the left panels as the red and blue lines respectively. 



the observed ridge profiles obtained from PICARD/SODISM and SOHO/MDI intensity data 
are compared to the simulated profiles. From these simulations, we are able to understand: 

• why the peak in the V-v diagram is shifted with / (the shift with the original l-v ridge 
increases as I increases); 

• why the peak gets closer to the original ridge as we get closer to the limb; 

• why the V-v ridge leaks towards higher frequencies. 

Figure [6] (left) shows that the energy in the p-mode range peaks at the extreme limb, as 
was already seen using MDI intensity images [15] or the guiding pixels of SOHO/LOI [2]. If 
we normalize this by a proxy of the noise level obtained by integrating the energy in the low- 
frequency band (right panel), we can see that we effectively have an amplification of the signal- 
to-noise ratio at the extreme limb for the three instruments. This peak is however slightly shifted 
inward for SODISM. This may be an instrumental effect but this needs further investigations. 

3. Conclusions 

The first helioseismic analyses of SODISM full- continuum images have been carried out for both 
the medium-/ and limb programs. We fitted mode parameters for a wide range of eigenmodes. 
The fitted peak asymmetries and frequency splittings are consistent with what is known from 
previous studies. We performed a first inversion of the rotational splittings and infered the 
internal rotation rate down to the radiative zone. The steep gradient of the tachocline is clearly 
seen at the base of the convection zone. The signature of p-modes can also be observed in limb 



1x10 7 
8x1 Cf 
| 6X10 6 

o 
z. 

a 

W 4x10 6 
2x10'"' 



[ 2700.21 - 
[ 2800.51 
[ 2900 82 
[3001 12- 
[3101.43 ■ 
[ 3201 .73 ■ 
[ 3302.04 ■ 
[ 3402.34 
[ 3502.65 
[ 3602 95 



2800.51 ]ti Hz 
2900.82 ] ft Hz 
3001 12] ft Hz 
3101,43 ]tt Hz 
3201.73 ]// Hz 
3302 04 J ft Hz 
3402.34 ] ft Hz 
3502.65 ] ft Hz 
3602.95 ] ft Hz 
3703 26 ] ft Hz 



i i i i i i i i i l I i i i i i i i i 




D.97 



0.98 



0.99 
R/Rlimb 



1.00 



2 0.4 



1.01 





i 1 > i 

— HMI ". 




MDI 


I y 1 1 


\ PICARD 


■ i < 1 1 < i i , , 





1.00 

r/rjimb 



Figure 6. Left: spectral energy level integrated in several frequency ranges of the 5-min 
oscillations envelope. Right: energy of the 5-min envelope normalized by the low frequency 
energy, for SODISM, SOHO/MDI and SDO/HMI. As far as SODISM is concerned, the 
amplification curve is displaced toward the inner part of the limb. This could be due to an 
instrumental effect not well taken into account in the data reduction pipeline. Note the narrower 
profile in the case of HMI, consequence of a better PSF. 



data and a signal amplification was detected at the extreme limb. By computing the response of 
any mode to the limb transform, we explain most of the features of the power spectra observed 
at the limb when compared to spectra built using full images. However, we have to face strong 
instrumental and orbital effects that affect both the photometric signal and the geometry of the 
solar images. We know from the past that our capacity in assessing the reliability of the details 
of many features analyzed using helioseismology (tachocline shape, polar jet, meridional flows, 
surface effects for structure inversions, etc. see e.g. [13J) greatly benefit from the cross-analysis 
made using different instruments (e.g. GONG and SOHO/MDI). We are therefore still working 
on improving our calibration procedure and we hope that PICARD will be extended at least 
one year in order to gather more intensity data that can be compared and cross-calibrated with 
SDO/HMI intensity and velocity images. 
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